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Abstract 

We determine the Tevatron's reach in supersymmetric parameter space in trilepton, like- 
sign dilepton, and dilepton plus tau-jet channels, taking all relevant backgrounds into 
account. We show results for the minimal supergravity model. With a standard set of 
cuts we find that the previously unaccounted for W^* background is larger than all other 
backgrounds combined. We include cuts on the dilepton invariant mass and the VF-boson 
transverse mass to reduce the W'j* background to a reasonable level. We optimize cuts 
at each point in supersymmetry parameter space in order to maximize signal-to-noise. 
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1 Introduction 



Low-energy supersymmetry (SUSY) is the most popular extension of the Standard Model 
(SM). It is vigorously sought at LEP, and has been and will continue to be actively looked for 
in the previous and forthcoming runs of the Fermilab Tevatron Collider 

There are many possible manifestations of low-energy SUSY. With more than 100 new 
parameters, theorists have out of necessity invented high scale models with drastically fewer 
parameters. These models can have qualitatively distinct low-energy spectra, leading to a 
variety of collider signatures. In this paper we explore the reach of the Tevatron in the 
parameter space of the most commonly considered model, the minimal supergravity model 0. 
As in many models, this model respects gaugino mass unification. This implies that in the 
physical low-energy spectrum the electroweak gauginos are significantly lighter than the gluino, 
so that their production cross sections are the largest in the allowed regions of parameter 
space. In addition to a large production cross section we want a significant branching fraction 
into a channel with relatively small Standard Model background. With all this in mind, the 
trilepton (3L) signal, l^l^l" Ipr with £ = e or /i, has been considered a gold-plated mode 
for SUSY discovery at the Tevatron ||^, |, ^ ^, 0, ^, |^, prompting several Run I analyses at 
the Tevatron [T^. The 3L signal is mainly produced via pp — > '^Xi ■ Being one of the most 
extensively studied channels for SUSY discovery, it was naturally among the main focal points 
of the Run II Workshop where the emphasis was placed on optimizing the analysis cuts in 
order to maximize the Run II Tevatron reach. 

In a recent paper p|, we took this approach further by considering thousands of sets of 
cuts in order to determine which one gives the best reach. We also supplemented our trilepton 
SUSY search with two other promising signatures — the inclusive like-sign dilepton |Tl| and 
'dilepton plus a tau jet' channels. Along with the mandatory plots of the Tevatron reach 
in parameter space, the main result from ^ was that the SM background has been grossly 
underestimated in the previous studies and to some extent in 0). We traced the main 

cause of the problem to the inadequacy of the event generator ISAJET [|T^ in simulating the 
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SM WZ and ZZ backgrounds. In ISAJET the zero width approximation is used in generating 



both WZ and ZZ. In PYTHIA [0, a Breit-Wigner distribution is used for the W- and Z- 
bosons. The finite Z-width leads to broader dilepton spectra and hence significantly larger 
background. 

The W'-f* background is not incorporated in either ISAJET or PYTHIA. Hence, it has not 
been taken into account in previous studies (see, however, Refs. [0, |^). We find that this 
background is larger than all previously considered backgrounds combined. In light of the 
importance of the trilepton channel for Run II, we are compelled to update our analysis of 
Ref. i. 

We were faced with several options as to how to incorporate the W'y* processQ There are 
several parton level Monte Carlo generators which use the full set of diagrams (see Fig. [^) to 
generate what is loosely called 'WZ", but in reality is the 2 — 4 process pp — > i^ui'^i''. 



Three such generators are MADGRAPH COMPHEP ^ and MCFM The choice of 
a particular generator is dictated by a matter of convenience and/or experience. We want to 
not only generate WZ events with the correct kinematics, but also to include a full detector 
simulation as we did in |^, making use of the SHW package [0, |TB|, |T^. In addition, to make 
the simulation fully realistic, we need to include the effects from initial and final state radiation 
(ISR,FSR), therefore we cannot just link one of the leading order parton level Monte Carlos 
to our detector simulation package. What we choose to do instead is to use COMPHEP to 
generate hard scattering events at leading order, then we pipe those through PYTHIA which 
adds showering and hadronization, and finally we run the result through SHW [|. The resulting 



parton-level cross section was integrated with the CTEQ4m structure functions [pT 



Unfortunately, with a standard set of cuts the WY background is about 2.7 fb, which 



*By Wj* we implicitly refer to the Z-j* interference as well. 

'''Alternatively, one can omit the first step and generate the WZ events directly from PYTHIA, reweighting 
the events so as to fit the distributions of a few key variables (e.g. dilepton invariant mass, lepton pT spectrum 
or angular distributions, etc.). In the early stages of this project we followed this approach and reweighted the 
PYTHIA events to fit the invariant mass distribution from MCFM. 'We then applied the same cut optimization 
procedure as in and presented our results for the Tevatron reach in a series of talks |po| . This procedure is, 
of course, only an attempt to approximate what we are doing here. The results turn out to be in reasonable 
agreement with the current results. 
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Figure 1: The diagrams for the pp — > VF+(Z/7*) i'^v^ii^l background. Here 
u and d stand for a generic up-type and down-type quark, respectively. 



is larger than all previously considered backgrounds combined (2.1 fb P]). This new source of 
background dwarfs previous estimates. For example, it is over 4 times the total background 
found in Ref. 0, and our total background is now more than 8 times the total background 
reported in 0. 

These recent developments necessitate the invention of new cuts, specifically designed to 
suppress the off-shell Z/7 component of the background. One obvious variable to consider 
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Figure 2: The invariant mass distribution rrii+i- of the opposite sign, same flavor 
leptons in (a) /i^e"'"e~ events and (b) e^e~^e~ events. The histograms show the 
results from COMPHEP and from PYTHIA (shaded). We have imposed nominal 
charged lepton cuts pt{^) > 5 GeV and m£+£- > 10 GeV. Each histogram is nor- 
malized to its cross section. In (b), we fill both invariant mass combinations, each 
with weight 1/2. 

is the invariant mass rrii+i- of an opposite sign, same flavor lepton pair in the event. The 
inclusion of the off-shell photon contribution increases the relative weight of events with low 
m£+£-. In anticipation of this effect, in Ref. we employed a low invariant mass cut of 
rrii+i- > 12 GeV. 

In Fig. H we show the mf+^- distribution in W Z events from COMPHEP and PYTHIA, 

before detector simulation and without ISR/FSR0. We divide the WZ trilepton sample into 

opposite flavor (OF) (e^/i"*"/!" and /i^e"'"e~) and same flavor (SF) {js^ jj,'^ ji" and e^e'^e~) 

subsets, and show the results for each subset separately in Fig. 0a and Fig. ^b, respectively. 

For the OF sample, we know unambiguously which two leptons came from the off-shell 2'/7, 

so we enter one invariant mass combination per event. However, in the SF sample, there are 
■fin ISAJET the invariant mass distribution of the leptons is a (5 function at the Z-mass. 
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two possible invariant mass combinations for each event, and there is no way to know which 
one was from the Zj^. Hence, in Fig. we enter both combinations, each with weight 1/2. 

First we see that neglecting the virtual photon contribution and the Z— 7 interference leads 
to a significant underestimate of the W Z background. In fact, the virtual photon contribution 
diverges in the limit m£+£- 0! Second, the low-end invariant mass cut mi+i^ > 12 GeV 
that we used in is clearly not very efficient in suppressing the additional 7* background and 
the cut threshold needs to be increased. The optimum threshold will depend on the signal 
distribution, whose shape is controlled by the value of the chargino mass m-+ and is thus 
parameter space dependent. We therefore incorporate the low-end invariant mass cut into our 
optimization scheme, and we consider the cuts 'mJ+^_ > {10-60} GeV, in 5 GeV increments. 
We choose the optimal one at each point in SUSY parameter space (for further details on our 
optimization procedure, see 0). 

In Fig. ^ we compare the COMPHEP and PYTHIA pt distributions of the leptons. We 
see that most of the additional events due to the 7* contribution tend to have small pt- This 
implies that the soft cuts on the lepton pt introduced in Ref. may be inefficient in removing 
the new background component. The soft pt cuts could be detrimental to the reach in regions 
of parameter space where the size of the background is important. 

Alternatively, Ref. suggests a cut on the transverse mass rriT of any iu pair which 
may originate from a IV-boson. The advantage of this cut is that it removes background 
events irrespective of whether the remaining lepton pair came from a Z, 7* or the interference 
contribution. We shall therefore optionally incorporate this cut in our analysis of all three 
channels: 60 < mT(£, z/) < 85 GeV. The remaining cuts that we use are fully described in 
Ref. [§ and will not be repeated here. 

We present our results for the Tevatron reach in the trilepton, like-sign dilepton and dilep- 
ton plus tau jet channels in Figs. ||, ^ and respectively. We require the observation of 
at least 5 signal events, and present our results as 3cr exclusion contours in the Mq — M1/2 
plane, for two representative values of tan (3, 5 and 35. We fix /i > and Aq = 0. The cross- 
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Figure 3: The same as Fig. ^, but for the px distribution of the leptons possibly 
coming from the Z. In the case of //^e"'"e~, we fill the pT of both and e~, each 
with weight 1/2. For the case of e^e^e" we fill the px of the odd-sign lepton with 
weight 1/2 and the px of the like sign leptons with weight 1/4 each. 

hatched region is excluded by current limits on the superpartner masses. The dot-dashed 
lines correspond to the projected LEP-II reach for the chargino and the lightest Higgs masses. 
In Figs, (a) the left dotted line shows where m,>^ = m-± and the right dotted line indicates 
TUf^ = '>TL^± (and rrif ~ ~ mg). In Figs, (b) the dotted lines show where mg^ = m-± (left) 
and = m-± (right). 

We see that although the inclusion of the 7* background leads to a significant increase 
in the raw background cross section, the reach is somewhat similar to what was presented in 
Ref. [Q, since the additional cuts help to increase the signal-to-noise ratio reasonably close 
to previous levels. At small tan f3 the trilepton channel provides for significant reach at both 
small Mq {Mq < 150 GeV) and large Mq {Mq > 400 GeV). The other channels have somewhat 
less reach. At large tan f3 the dilepton + r jet channel provides the best reach at small Mq 
(Mq < 160 GeV), while at large Mq {Mq > 400 GeV) the trilepton channel still provides for 
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Figure 4: Tevatron reach in the trilepton channel in the Mq — plane, for fixed 
values oi Aq = 0, fj, > and (a) tan/3 = 5, or (b) tan /3 = 35. Results are shown for 
2, 10 and 30 fb~^ total integrated luminosity. 



decent reach with 30 fb~^. With only 2 fb~^ the reach is quite hmited. 

In Fig. 1^ d^) we show the optimum cuts chosen in our optimization procedure, in the Mq, 
Ml/2 plane, for tan/? = 5 (tan/? = 35), in the small Mq region. We use the following notation 
to describe the set of cuts at each point. The central symbol indicates the set of lepton px 
cuts: the symbols "1" through "5" refer to {11,5,5}, {11,7,5}, {11,7,7}, {11,11,11} and 
{20, 15, 10} GeV lepton pt cuts, respectively. The left superscript shows the value (in GeV) 
of the low-end invariant mass cut (mJ+^_ > 10 to 60 GeV). A left subscript "T" indicates 
that the cut on the transverse iu mass was selected. The right superscript shows the 
cut: > {15,20,25} GeV ("15", "20", "25"), or no cut (no symbol). A right subscript 
denotes the high-end dilepton invariant mass cut: \m£+£- — Mz\ > {10, 15} GeV ("10", "15") 
or me+e- < {50,60,70,80} GeV ("50", "60", "70", "80"). And finally, a tilde over the central 
symbol indicates that the luminosity limit came from requiring 5 signal events rather than 3a 
exclusion. 
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Figure 5: The same as Fig. ^ but for the like-sign dilepton channel. 



In Fig. 1^ we see that in the regions where background is an issue, the combination of the 
mT cut and a tighter low-end dilepton mass cut TnJ+^_ ~ 20 GeV is typically preferred. Indeed, 
we find that these additional cuts reduce the WZ background by more than a factor of 3, from 
4.1 fb (with soft cuts 0) to 1.2 fb. Notice, however, in the small tan/? case the transverse 
mass cut is never enough by itself, i.e. whenever it is chosen, it is almost always supplemented 
with a rnj+^- cut of 15 to 25 GeV (with the exception of two points with high lepton px cuts). 
On the other hand, there are significant regions where the low invariant mass cut rnj+^^ by 
itself is enough to kill the background, and the transverse mass cut is not needed. In the large 
tan/3 case the transverse mass is always chosen at small Mq, but only occasionally at large 
Mo. 

We should point out that the optimum cuts in Figs. ^ and || can be interpreted in two 
ways. First, for a given total integrated luminosity, say 10 fb~^, one can first roughly look up 
from Fig. | the sensitivity reach of the Tevatron. Then, for the parameter space well inside 
the sensitivity region, the actual choice of cuts is not so crucial. However, as one approaches 
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Figure 6: The same as Fig. but for the dilepton plus a tau jet channel. 



the boundary of the sensitivity region, the choice of optimum cuts as a function of parameter 
space (as opposed to a fixed, non-optimized set of cuts) can enhance the reach by an additional 
10-20 GeV along the M1/2 direction ||^. Alternatively, at a given parameter space point near 
the border, optimization can reduce the total integrated luminosity required to observe or 
exclude that point by up to a factor of two 0. 

In conclusion, we find that the trilepton channel remains one of the leading candidates 
for SUSY discovery at the Tevatron. The other two channels are in a sense complementary, 
although not as powerful. The dilepton plus tau jet channel can be combined straightforwardly 
with the trilepton channel to maximally increase the reach. With the new very important 
PF7* background included, the reach of course suffers somewhat. We find that with only the 
standard Run 2 luminosity of 2 fb~^ the reach is quite limited. With the larger background 
it is even more imperative that the Tevatron collect as much luminosity as possible to have a 
decent chance at discovering supersymmetry. 
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Figure 7: The optimal sets of trilepton cuts in the Mq, M1/2 plane, for tan/3 = 5 
and small Mq. We show the optimal low end dilepton mass cut mJ+^_, missing Et 
cut high end dilepton mass cut rrii+i-, transverse mass cut and lepton 
cut (see text). The dotted lines indicate the reach contours from Fig. ^. 
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Note added: Our results shown in Figs. are similar to those of Ref. once the error 
in the original version of ^ was fixed We warn the reader that the analysis of still 
neglects the 7* contribution and the interference effects in ZZ production, and lacks a detailed 
detector simulation for the W Z process. 
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